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I. INTRODUCTION 



Ultrarelativistic heavy-ions experiments at CERN and Brookhaven have provided im- 
portant information on the processes involved in high-energy nucleus-nucleus collisions. 
Global features such as stopping, multiplicity, transverse energy, and rapidity distribu- 
tions of particles can be described to a good approximation by models based on multiple 
hadronic collisions. However, several results in central nucleus-nucleus collisions, such 
as strangeness enhancement J enhancement ||, and large fluctuations in multi- 
plicity and transverse-energy distributions 0|| , deviate from naive extrapolations based 
on summing nucleon-nucleon or nucleon-nucleus collisions. The crucial issue is to deter- 
mine whether these interesting signals arise from the presence of a quark-gluon plasma 
or from more conventional phenomena such as rescattering and thermalization processes 
in the hot hadronic gas formed in the wake of the collisions. As we describe here, study 
of fluctuations provides a particularly valuable probe of the physics underlying these 
processes. 

An early indication of the importance of understanding fluctuations in observables 
was provided by the NA34 geometrical model This model gave reasonably accurate 
parametrizations of the NA34 multiplicity and transverse energy (Et) spectra in terms 
of summing multiple nucleon-nucleon collisions - each of which provides a distribution of 
transverse energy and multiplicity - with the collision geometry incorporated. However, 
as Baym et al. || pointed out, the fluctuations extracted in this way from experiment, 
as reflected in the widths of the high-energy tails of Et distributions, are much larger 
than those expected from a microscopic model with the known fluctuations in transverse 
energy in nucleon-nucleon collisions. 

We focus here on understanding the subtle and interlocking roles of the nucleon cor- 
relations in the underlying nuclear structure of the colliding nuclei, the range of the 
high energy nucleon-nucleon interactions, and cross-section fluctuations in determining 
fluctuations in observables. In Refs. and ]7| we gave a brief account of the importance 
of including cross-section fluctuations and nuclear correlations in calculating Et spec- 
tra. Calculations with these effects included showed that one should expect the large 
fluctuations apparent in the data. The present work includes a more detailed derivation 
of these effects. 

Short-range correlations, which arise mainly from the short-range repulsion between 
nucleons as well as Pauli blocking, exclude other nucleons within a distance of order 
the correlation length, and reduce the fluctuations around the mean number, (N), of 



zone is small compared to the volume of the nucleus but large compared to that of a 
single nucleon. For a target nucleus containing an uncorrelated gas of nucleons freely 
running about, the number of target nucleons in the interaction zone can range from 
none to all, as given by the binomial distribution, leading to large flucutations in N. 
However, the nuclear matter in the nucleus is in fact a strongly correlated self-bound 
liquid, not a gas; the density fluctuations are shorter ranged, and the fluctuations in the 
number of nucleons in the interaction zone greatly reduced. This effect is even more 
significant for nucleus-nucleus collisions. 

Since the correlation length is of order the internucleon distance, itself of order the 
range of the interactions or the cross-sectional distance r s = (cr/V) 1 / 2 , fluctuations in 
observables depend as well on the detailed interplay of the correlations and the range 
of the interactions. A normal assumption in microscopic models is that the projectile 
interacts with all nucleons in the tube along its trajectory with area given by the total 
NN cross section. However, allowing a more diffuse interaction range tends to reduce the 
importance of correlations, and increase the fluctuations, as will become evident below. 

We primarily develop the study of correlations and fluctuations within the frame- 
work of binary collisions of target and projectile nucleons, but we briefly look as well at 
the predictions of the "wounded-nucleon" description of sources in terms of participant 
nucleons. The underlying basis of the binary-collision picture is the cutting rules in the 
field-theoretic diagrammatic analysis of high-energy scattering processes that were found 
by Abramovskii, Gribov, and Kancheli (AGK) || and Bertocchi and Treleani Q within 
the inelastic eikonal approximation. These rules imply exact cancellations between in- 
elastic interactions of the projectile with more than one nucleon at a time. They thus 
imply that in nucleon-nucleus scattering, away from the projectile fragmentation region, 
inclusive particle production is that in NN scattering times the number of nucleons, A, 
in the target nucleus; in nucleus-nucleus (BA) scattering away from the nuclear fragmen- 
tation regions, the result is that in NN scattering times BA, the numbers of nucleons in 
the projectile and target. In the fragmentation regions, shadowing effects must also be 
taken into account. In the Appendix we discuss the kinematical regions where the AGK 
theorem is valid; we also indicate how one can establish the correspondence between the 
basic probabilistic formulae deduced in this paper from geometrical arguments and the 
field-theoretic combinatorics of the AGK approach.! 



2 The experimental situation does not make a clear case for either an underlying binary- 



Cross-section fluctuations also increase the fluctuations in the number of collisions, 
as shown in ||. The basic picture is that, owing to frozen color configurations in the 
projectile, the probability of a given cross section, a, in a primary collision is given 
by a distribution function, -P(cr), instead of the cross section always taking a constant 



mean) value. In Ref. [0 we discussed in detail the nature of cross-section fluctuations in 



nucleon-nucleon interactions, and for pion-nucleon interactions in Ref. ||13|| , and showed 
how the distribution function P for cross section fluctuations can be determined from 
both inelastic shadowing in nucleon-deuteron scattering, and from diffractive excitation 
in high-energy nucleon-nucleon scattering. In the Appendix we briefly discuss cross- 
section fluctuations within the AGK framework. 

This paper is organized as follows. In Sec. 2 we first calculate the effect of nuclear 
correlations on fluctuations in the number of multiple NN interactions in proton-nucleus 
collisions. We then generalize the results to nucleus- nucleus collisions in Sec. 3. In the 
following section, 4, we extend the analysis to include contributions from cross-section 
fluctuations, and then, in Sec. 5 we apply these results to the question of the fluctuations 
in the tails of multiplicity and transverse energy distributions. In the final section we 
summarize and outline the procedure necessary to implement the effects discussed here 
in the Monte-Carlo descriptions of AA collisions, and in the Appendix discuss the AGK 
theorem. 



ing) as well as coherent processes, e.g., showers of particles overlapping spatially and interact- 
ing strongly, which the binary picture does not account for. Proportionality of the inclusive 
spectrum to A at mid-rapidities is observed experimentally, e.g., in the recent systematic 



study by Whitmore et al. of pA scattering at 100 and 320 GeV lab energy [1C|. The faster 
than ~ A dependence for low total momentum secondary hadrons observed near the nuclear 
fragmentation regime is likely due to secondary hadron interactions, including interaction of 
produced hadrons with residual nucleons. On the other hand, experimental 16 0+Em(ulsion) 
data [|ll]] indicate that total multiplicity is intermediate that of the binary-collision model and 
the wounded-nucleon model. 



II. CORRELATIONS AND FLUCTUATIONS IN PROTON-NUCLEUS 



COLLISIONS 

To gain a first understanding of the effect of nuclear correlations on fluctuations we 
begin by examining collisions of protons with heavy nuclei, and treat the more compli- 
cated case of nucleus-nucleus collisions below. Neglecting cross-section fluctuations at 
this point, we take the NN inelastic cross-section to be its average value, a = a ~ 32mb 
(at CERN energies). 

The mean number of nucleon-nucleon collisions, N, that a proton at impact parame- 
ter b makes with the target nucleus can be written, in the impact parameter formalism 



of Ref. 0, as 

(N) = Jd 3 rp A (r)^- b (b-r ± ), (1) 

where pA is the nuclear density distribution, r_|_ is the transverse coordinate, and 
(da/d 2 b)(b — r±) is the probability for the projectile nucleon with impact parameter 
b to collide with a target nucleon at position r. In the eikonal approximation! with 
straight-line geometry for multiple NN collisions - the Glauber approximation - da j d 2 b 
depends only on the transverse coordinates r_|_; it has an interaction range ~ r s ~ 1 fm 
and obeys the sum rule 

/A-4(rJ = „. (2) 

As we shall see below, it is important to take the range of the interaction into account 
in the multiple NN scattering problem since it is comparable to the average internucleon 
distance r = (3/47rpo) 1 ^ 3 — 116 fm (where po — 0.16 fm -3 is nuclear matter density), as 
well as the correlation length, £, between nucleons in a nucleus, ~ 1 fm. Since r s <C Ra, 
the target radius, the mean number of NN interactions can, aside from a small correction 
at the nuclear surface, be calculated without detailed knowledge of da/d 2 b; using (0) 
we have 

/oo 
dzp A (b,z)=oT A (b), (3) 
-oo 

where T A (b) is the thickness function for the target nucleus with mass number A. If 
the nucleus is assumed to have constant density p with radius R A = r A 1//3 , then 
T A (b) = 2p (R A - b 2 f 2 . 



The fluctuation in the number of collisions, N, at impact parameter b, is similarly 
given byi 

(N(N - 1)) = / d s rd 3 r'pf(T,r')^- b (b - rj^(b - r' ± ); (4) 

(2) 

here p A is the two-body density distribution 

p? ) (r,r / ) = ^ t (r)^ t (r , )^(r , )^(r)) 

= p A (r)p A (r')(l-C A (vy)); (5) 

where if)^ and if) are the nucleon creation and annihilation operators. Because nucleons 
tend to stay away from each other at short distances, each nucleon is surrounded by a 
correlation hole, described by the correlation function C A (r, r'), which generally, as a 
function of r', is non- vanishing about the nucleon at r over a range £ in |r — r'| of order 
1 fm in nuclear matter. From the commutation relations of if) and if)* one readily derives 
the sum rule on the correlation hole 

Jd 3 r'p A (r')C A (r,r') = l, (6) 

which says that the hole depletes precisely one nucleon. Equation fl5|) implies that 
C A < 1. For a nuclear liquid, one expects C A to be non- negative, as is borne out by 
detailed calculations |15|.i 

With (Jl), we can write (f|) as 

(N(N - 1)) = J d 3 rd 3 r'p A (r)p A (r')(l - C A (r, r'))^(b - r^)^(b - r' ± ) 

\N) 2 - (N)a) , (7) 

where we define the dimensionless parameter a by 



4 In writing (|j) we neglect Glauber shadowing in the interaction of the projectile with several 
nucleons of the nucleus. As discussed in the Appendix, AGK cancellation H|,[| in the Glauber 
model justifies this neglect when calculating average characteristics of the final state away from 
the projectile fragmentation region. 

5 One can imagine situations in which C A could be negative, e.g., halo nuclei ( 11 Li) or helium 
clustering inside light nuclei. In these systems high density "lumps" inside the nucleus increase 



(N)- 1 J d 3 r'p A (r)p A (r')C A (r; r')^(b - rj_)^(b - r' ± ). (8) 



a 



Physically, the parameter a measures the fraction of the overlap between the collision 
volume and the volume of the correlation hole described by C A (Fig. 1). To see this let 
us take the correlation function to be spherically symmetric, i.e., C A (r, r') = C A (\r — r'\). 
Noting that for large nuclei both r s and £ are much smaller than the nuclear radius, and 
for non-peripheral collisions the integral in (|j) is dominated by the bulk of the nucleus 
and is independent of impact parameter, we see that Eq. (||) reduces to 

do , N da 



a 



& J dzd*r ± cPr> x C A ( V^ 2 + (r ± - ^f) ^(r ± )^(rl) . (9) 



From Eq. (|7|) we obtain the variance in the number of collisions, 

UwS <5M = 1 _ a . (10) 

This equation shows that nuclear correlations, through the quantity a, reduce the fluc- 
tuations in proton-nucleus collisions; we shall see below that nucleus-nucleus collisions 
are even more sensitive to nuclear correlations. The parameter a varies from zero, for an 
uncorrelated gas of nucleons, to unity, for a strongly correlated system, independent of 
the detailed forms of the correlation and collision probability functions. It is instructive 
first to evaluate a in various limits, and then consider the general case. 



A. Extreme optical limit 

In the extreme optical limit in which the range of the interation r s is negligible, i.e., 
r s <C £, the NN scattering probability, da/d 2 b, can be approximated by 

^(b-r ± )^ 2 )(b-r ± ). (11) 
d z b 

Then from Eq. @ we find that 

a ~ ap J dzC A (z) , r s < £. (12) 

The parameter a can thus be interpreted as the optical depth of the correlation hole, the 
width in nucleon mean-free paths, l/o~p A , across the hole. 

If the correlation function is a simple step function, C A (\r — r'|) = 0(£ — |r — r'|), 
then the sum rule (§) implies that £ = tq = (3/47rpo) 1//3 ; and we obtain 



B. Ideal gas 



In an ideal gas with A distinguishable particles, the correlation function is a constant, 
Ca = 1/A, determined by the sum rule; the correlation "hole" extends without structure 
across the entire nucleus (see Fig. 2). In this case we can also use Eq. ([El]) and find 
from Eq. (H) that for collisions at small impact parameter, 

a=(N)/A, (14) 

which coincides with the result ( |i~3l) for £ = R. The term 1-a in Eq. ( |10D is then the 
standard binomial result for the fluctuations in the number of collisions of a projectile 
passing through an uncorrelated nuclear gas, with individual collision probability p = 
(N)/A = a. This is the minimum value for a, corresponding to maximum fluctuations, 
and correlation length £ of order the nuclear radius. 



C. Maximum correlations 

When ^ < r s we can write 

p A {r)C A (r,r')~5®(r-r'), (15) 
and find from Eq. @ that 

a= \! d2h {ik {h) ) ' e<<rv (16) 

If we assume, as one often does in event generators simulating multiple NN collisions, 
that the projectile nucleon collides with all nucleons situated within a tube of cross- 
sectional area a along the beam axis with impact parameter b (see Fig. 1), i.e., 

^(b-r ± ) = 6(r s -|b-rj), (17) 

we find 

a = l, £<r s , (18) 

the maximum value of a. As Eq. fllTf) shows, the fluctuations in the number of binary 
collisions vanishes in this case. This result is due to the fact that in this limit each 



D. Gas with excluded volume 



In event generators, nucleon correlations rule approximated by not allowing 

two nucleons to be situated within a distance, £, of each other. The correlation function 
becomes in this way that of a gas with an excluded volume (see Fig. 2), which to a first 
approximation, when £ < r , is given by: Ca = 1 for r < £ and, from the sum rule, 
Ca — ( r o — £ 3 )/Ra f° r r > £• Since a is linearly proportional to the correlation function, 
it is the weighted sum of the two pieces. To evaluate the first term, we use Eq. ([18]) 
when £ <C r s , and Eq. (|1^) when £ S> r s ; for the second term we use Eq. (|T4| ) generally, 
and find in the two extreme limits, 

1 i»r,. 

A numerical calculation for general £/r s with the "tube" cross-sectional area of Eq. fll7|) 
results in fluctuations well-approximated by the first line of Eq. ([I9|), except for larger 
values of £, where a is reduced by approximately 20%. It is important to note that 
the fluctuations and a are very sensitive to the choice of £. In the following subsection 
we treat more generally the situation where the range of the correlation function is 
comparable to r s . 



E. General case 

The nucleus is a saturated self-bound liquid and is thus strongly correlated, as can be 
seen from the correlation function of Fig. 2 taken from Ref . [15j . Since the correlation 
length £ is of order the internucleon distance r = (3/47rpo) 1 ^ 3 — 1.17 fm, itself of order 
r s ~ 1.0 fm, a depends on the detailed structure of Ca and da/d 2 b. To make the 
discussion of the general situation more specific, we construct simple parametrizations 
of both the correlation function and the scattering probability, da/d 2 b. In particular we 
assume a Gaussian form for the correlation function (see Fig. 2), 

C, (r y) = e-(r-r<)72? 2 ; (19 ) 

the sum rule (|6[) relates the correlation length £ to the nuclear density by 

£ = r (2/97r) 1/6 ~ 0.75 fm . (20) 



In the impact parameter formalism of Ref . |14| , the imaginary part of the scattering 
amplitude, Im f(q), for momentum transfer q, is related to the collision probability 
function by 

W(g) = ^/A- b ^(b), (21) 



where k is the incident momentum. Note that as q — > 0, Eq. fl2T|) yields the optical 
theorem, 

a = yIm/(? = 0) . (22) 
At high energies, the imaginary part of the scattering amplitude can be parametrized as 

Im/(g) = ^e-«/ 4 ; (23) 
the slope parameter r 2 depends weakly on energy, and at CERN energies, Ei a b ~ 200 



GeV, r a ~ 1.2 fm ||14|| . With this form we find 

da „ . a _ft2 / r 2 

— b =— e ^ 24 

Substituting the two Gaussian parametrizations, ( |T§| ) and ( P4"| ) in the definition of a, 
Eq. (W), and making use of the fact that both £ and r s are i?, we find 



Po<? 

a = 



J ^e- 2 / 2 « 2 |d 2 r ± d 2 rle-^- r 'x) 2 /^ e -(^ 2 +^ 2 )/^ ; (25 ) 



which, by transforming to coordinates ± r^_, and using the relation ( P0{ ) is readily 
evaluated as 

a = ieh- (26) 

In order that the collision probability function, (§J), be less than unity, the condition 
a < 7rr 2 or equivalently r s < r a must be satisfied. (Note that r s = (a/ir) 1 / 2 determines 
the strength of the interaction whereas r a determines its range.) For the values for r a 
and a quoted above a < 7rr 2 is indeed satisfied. This condition, used in fl2"6f) , implies 
that a < 1/2. This surprising result does not in fact disagree with the earlier result 
a = 1 in (|18D since there da/d 2 b was assumed to be a step function. Employing a 
diffuse form for da/d 2 b has the physical consequence that two nucleons lying on the 
same projectile trajectory are not necessarily both hit, which softens the restrictions 



F. Numerical studies 



Pieper [[U| has calculated the detailed probability distribution Pat(6) for a proton 
colliding with various closed-shell nuclei numerically using correlated many-body wave 
functions, with the assumption that da/d 2 b is given by the step function of Eq. (|17D. The 
distribution function is considerably narrower for central collisions than a Poisson. By 
counting the number of nucleons in cylinders of radius 1 fm (corresponding to o = 32 
mb) he finds, for central events, a ~ 0.31 — 0.33 for both small (A=40) and large 
systems (A=224). In addition, the fluctuations are found to increase towards unity 
(corresponding to vanishing a) for peripheral events. [Taking only the correlations of a 
free Fermi gas of neutrons and protons into account yields a = 0.25 for central collisions.] 
Pieper's value for a is larger than found above, a = 0.25, primarily because of his use 
of a step-function form for da/d 2 b instead of a more diffuse Gaussian. 



III. FLUCTUATIONS IN NUCLEUS-NUCLEUS COLLISIONS 

We now generalize the calcuation of the number of NN collisions and its fluctuations 
to nucleus-nucleus collisions. We continue to assume the Glauber approximation with 
straight line geometry; the assumption that the nucleons do not move between successive 
hits is valid at high energies as a consequence of Lorentz contraction of the nuclei and 
the relatively small transverse momentum transfers in NN collisions. 

The mean number of NN collisions in a central (b = 0) collision of a projectile nucleus 
B with a target nucleus A is given by 

(N) = / d 3 rd 3 r'p B (r)p A (r')^(r ± - r' ± ) 
~ a j ' dz dz'd 2 spB(s, z)pa(s, z') 

= aj d 2 sT B (s)T A (s), (27) 

where we assume a <C Rb 2 ,R a 2 - When the projectile is small, i.e., Rb 2 <C Ra 2 , the 
mean number of collisions becomes 

(N) ~ Bap R A , B < A , (28) 

which is just B times the number of NN collisions in a central proton-nucleus collision. 
The fluctuations in the number of NN collisions are found from 



+ I d 3 rd 3 r'd 3 r"p {2 \v, r'K(0^;(^ - r±)^(*'± ~ r'l) 

+ J ^VrfVW>g ) (r,rOp! 4 2) (r'',r w )^(^-r'_[)^(r' ± -rT). (29) 

The four terms in Eq. ([29]) correspond to: 1) the same projectile nucleon colliding with 
the same target nucleon, 2) the same projectile at position r colliding with two different 
nucleons at positions r" and r'" in the target, 3) a nucleon at r" in the target being 
hit by the two different nucleons at r and r' in the projectile nucleus, and finally, 4) 
two different nucleons at positions r and r' in the projectile nucleus colliding with two 
different nucleons at positions r" and r'" in the target nucleus. 

The second term in ( p9|) is evaluated analogously to the proton-nucleus case and is 

= aj d 2 sT B (s)T 2 (s) - (N)a, (30) 

while the third term in ( |29] ) contains the two-body density distribution of nucleus B and 
gives a contribution 

= aj d 2 sT 2 B {s)T A {s) - (N)/3, (31) 

where (3 ~ a is the correlation parameter for nucleus B, defined analogously to Eq. (§). 
The fourth term in Eq. ( p9|) reduces to 



= (iV> 2 - a J d 2 sT B (s)T 2 (s) -a J d 2 sT 2 (s)T A (s) + >y(N) . (32) 

The four terms here arise from each nuclear two-body density distribution having two 
terms [cf. (H)]; they correspond respectively to the product: i) without any nuclear 
correlation functions, ii) with one correlation function from nucleus A, iii) with one 
correlation function for nucleus B, and finally iv) the product containing both correlation 
functions. The first three terms in (|32|), which contain only one or no nuclear correlation 
function, are straightforward to calculate and do not depend on the detailed forms of 
the correlation functions or the collision probability function, as long as £ and r a are 
much smaller than the nuclear radii. The last term is more tedious to evaluate because 
it contains the products of both correlation functions and two collision probabilities; 
these enter in the parameter, 7, defined by: 

7 = (N)' 1 f d 3 rd 3 r'd 3 r"d 3 r"' p B (r)p B (r')p A (r")p A (r'") 



xC B (r - r')C A (r" - r'")^± ~ r T)^( r ± " r T) ■ (33) 



The variance in the number of binary collisions, found from Eqs. (29)-(32), is 



uj n = 1 + 7 — a — (5. (35) 

This result generalizes the proton-nucleus result, Eq. (|10D , to nucleus-nucleus collisions, 
and improves the result derived in |J which assumed 7=1. With £ = 0.75, r a = 1.2 
fm and a = 32 mb we find a = (3 = 0.25, 7 = 0.15, and ujn = 0.65. 

The derivation of Eq. ( |35| ) assumes that the correlation length, £, is much smaller 
than the nuclear radii, Rb and Ra, as is the case for a strongly correlated nuclear fluid.! 
The importance of including the correct correlations is well-illustrated by calculating the 
fluctuations in the opposite limit in which the nucleons are allowed to swarm around in 
the nucleus as a free gas. The correlation functions are then Ca = 1/A and Cb = 1/-B, 
and when r a -C Rb <C Ra, one obtains directly from (f29"|) that 



(N 2 ) = (N) + a 2 Jd*sT B (s)n(s) (l - I 



1 \ / 1 \ / 1 

2 / j2„7i2/„\rr /„\ 1 1 \ , /Art 2 ' 1 



+a 2 y d'sn{s)T A {s) (l - - j + (iV)^ ^1 - - j ^1 - - j , (36) 

where the terms correspond to the successive terms in (p9"D- The resulting expression for 
the variance follows from Eqs. ( p9[) and (|36|): 



where we define 



c^at = 1 - H ~B~ NpB ( nuclear gas )' ( 37 ) 



N pA = (N)~ V 2 | d 2 S T B ( S )T 2 ( S ), (38) 
N pB = (iV)-V 2 / rf 2 sT 2 ( S )T A ( S ). (39) 



The quantities N p a and N p b are straightforward to calculate for small projectiles, 
since then Ta(s) ~ 2pqRa', N p a reduces to the number of NN collisions in central proton- 
nucleus collisions, 

N pA ~ 2ap R A , (40) 



and 



NpB ^ %jpqR b , B<^A. (41) 

The factor 3/2 instead of 2 follows from the impact parameter averaging necessary for 
the projectile nucleus in the limit B<A 

To calculate N pA and N p b more generally, we first assume uniform density nuclei 
with sharp surfaces. Then from ([]]) we find, 



(iV> = Anapl J" B db 2 (r 2 b - b 2 )'' 2 (r\ - b 2 ^ ' 



3 

-Bap R A 



R% db 2 [R 2 B - 

2 1 / „2 i\2i / £ + 1 



1 + x 2 (x 2 - lrln 



(42) 



2\ 



2x \x — 1 

with x = Ra/Rb- From (|38|) we have 

N pA =(N)- 1 a 2 J d 2 sT B (s)T 2 (s) = 8(N)- l a 2 npl f* db 2 (R 2 B -b 2 f /2 (R 2 A -b 

= A{N)-WBplR 2 A (l - 2 -x 2 ) , (43) 

and from fl39|), 

N pB = (iV)-V 2 / d 2 sT 2 B (s)T A (s) = 8(iV)-V 2 7r^ jf^ rf6 2 (i?| - 6 2 ) (i$ - b 2 f 2 

= A(N)- l a 2 BplR 2 B (x 3 - jj(x 5 - (x 2 - l) 5 / 2 )) . (44) 

For small targets, the assumption Rb <C -Ra used in deriving (|4"0|) and (|4T| ) overestimates 
and Ap£, as well as the fluctuations. For B = A, and b = 0, the above results give 

g 

N pA = N pB = -crp R A ■ (45) 

A diffuse nuclear surface decreases the average density, as well as decreases N pA and 
N p b from the values calculated above. A numerical calculation with a Woods-Saxon 
density distribution p = p /(l + exp (r — R)/S), where R is determined from the total 
number of nucleons in the nucleus and the surface diffuseness is chosen as 5 = 0.55, 
yields a further reduction of N pA and N pB of order 20% for A=2% to 10% for heavy 
nuclei, compared to the above results for a sharp surface. The effects are largest for 
small nuclei, where the nuclear surface is relatively large. 

The nuclear gas result (|3~7D is considerably larger than that of the strongly correlated 
nuclear liquid, Eq. (^), as can be seen in Fig. 3. The discrepancy can be traced 
back to the differences in the correlation functions. The derivation of (1351) assumed that 



P, e.g., as the value a = (N)/A in the proton-nucleus case, ([14]). The larger ojjq in the 
ideal gas case arises mainly from the second term in Eq. (|36|). The extra fluctuations are 
caused by "geometric correlations," that is, if one of the projectile nucleons collides with 
very few, or with numerous, target nucleons, then there is high probability that the other 
projectile nucleons will do the same since the projectile nucleons are nearby spatially. 
Correlations in nuclei, however, reduce the fluctuations in the number of proton-nucleus 
collisions and effectively cancel the geometric correlations. 

This exercise clearly shows the necessity treating nuclear correlations carefully in 
numerical simulations of relativistic heavy-ion collisions. In simulations the nucleons 
are often positioned successively in the nuclei by Monte-Carlo methods with the restric- 
tion that they should be at least a distance, £, apart. Such a procedure, as discussed 
above, describes a free gas with excluded volume, and fails to generate the correlations 
characteristic of a self-bound liquid. Generally, the correlation function resulting from 
this procedure has a correlation hole of radius £, but has a tail which oscillates for large 
distances and can be negative. While the correlation length, the range of the correlation 
hole, is close to the interparticle spacing in nuclear matter, £ ~ r , the minimum sep- 
aration distance of the centers cannot be directly translated into a correlation length. 
When £ > r , the integral of the density over the excluded volume exceeds unity, and 
the sum rule @ requires the correlation function to be negative on average outside the 
excluded volume. (Note that the maximum value of £ is ^ 2tq in order to accomodate 
all the nucleons in the nucleus). When £ < r the correlation function has a long tail 
characteristic of an interacting nuclear gas similar to the "excluded volume" correlation 
function of Fig. 2. 



For example, the Fritiof model requires the centers to be at least 1.13 fm apart fI8 



Numerically, this results in a ~ 0.4, and so for proton-nucleus collisions the fluctuations 
agree with those for the more-correct correlation function. However, in nucleus-nucleus 
collisions this failure to include liquid-like correlations generally overestimates the fluctu- 
ations, as Eq. ( |5TD illustrates. In the Fritiof model one finds, for example, that ojjq ~ 1.77 
in the case of 16 0+Au []rg| ], a value between that of a free gas, Eq. Q3"7|), and that for a 
nuclear liquid, Eq. (|3q) . The fact that a simulation reproduces large fluctuations does 
not imply that it correctly includes the physics producing the fluctuations. 

The density distribution resulting from positioning the nucleons by Monte Carlo with 
a minimum separation can also have the unpleasant feature that it is largest near the 
nuclear surface. This is because the imposition of a minimum distance between nucleons 
makes it harder to fit subsequent nucleons into the nucleus; it is easier to accommodate 



correct density distribution, but also describe nucleon-nucleon correlations correctly, 
guaranteeing that the correlation function is generally short ranged, vanishing within 
a distance less than the nuclear radius, and does not become negative. To implement 
correlations correctly in numerically Monte-Carlo codes requires constructing the initial 
nuclei by distributing all the particles in a nucleus at a time, rather than attempting 
to simulate correlations according to a rule that puts the particles in sequentially. The 
correct probability distribution for the nucleon positions is |\I/a(^i,^2, •••?"a)| 2 , where 
is the /L-particle ground state nuclear wave function. In practice, for given I^aP, one 
can follow a standard Metropolis algorithm to generate most-probable configurations. 
Many-body wave functions of Jastrow form, represented as the product of the mean-field 
wave function and two-nucleon correlators, readily lend themselves to development of 
generators of initial conditions for collision simulations.B 



IV. COLOR OPACITY AND TRANSPARENCY, AND CROSS-SECTION 

FLUCTUATIONS 

In hadron-nucleus collisions at high energies, the internal configuration of the color- 
carrying degrees of freedom of the projectile hadron are frozen by Lorentz time dilation 
fl9| , |20|| . The hadron- nucleon cross section is, as a consequence, characterized by a prob- 
ability distribution P(cr). The condition for a projectile nucleon to be frozen in passing 
through a target nucleus is that the crossing time, 2Ra/c, be less than the internal 



dynamic time, given by the uncertainty relation between time and energy |21 

2R A ^ = -^-^ = (46) 
lS.I1/ m N * — m N m^v* — ttln 

where is the nucleon mass, N* is the lowest-lying resonance with the quantum 
numbers of the nucleon, of mass m^* ~1.5 GeV, and 7 is the Lorentz factor. This 
condition is well-satisfied for projectile energies above p/ a f>~ 40 GeV/c for heavy nuclei, 



7 In particular, S. Pieper (private communication) of Argonne National Laboratory has gen- 
erated useful simplified approximations to many-body nuclear wave functions ( and un- 
published) of the Jastrow form that correctly reproduce the A-particle spatial distribution 
functions, although not the details of the spin correlations. These wavefunctions have recently 



as at the CERN SPS and LHC, and RHIC. The concept of cross-section fluctuations is 
well established in inelastic shadowing and diffractive hadron-hadron scattering, as well 
as in coherent diffractive-dissociation of hadrons scattering on nuclei, and is an essential 
feature of the dynamics of ultrarelativistic heavy-ion collisions fl6l,|T2 . 



The effects of cross-section fluctuations are several-fold. When a hadron is in a 
"small-sized configuration" its interactions are suppressed because of the small spatial 
extent of the color fields in the hadron, which leads to the phenomenon of color trans- 
parency |[22|| . On the other hand, a hadron in a frozen "large hadronic configuration" 
experiences a stronger interaction characterized by a cross section a > a, its average 
value, when passing through a nucleus, a "color opacity" effect. We now discuss, in 
further detail than in Ref. ||, and also for the more general case where the projectile B 
is not necessarily much smaller than the target A, how these effects lead to enhancement 
of the fluctuations in observables such as transverse energy and multiplicity. We then 
compare with fluctuations found at CERN. 

In calculating average quantities such as the mean number of collisions, (N), the 
cross section a can be replaced by its average value, a = (a)i ~ 32mb, where the 
subscript / refers to the average over internal configurations. However, in calculating 
fluctuations, one must first calculate for a given internal configuration and then average 
the final result. In particular the average of the square of the cross section for collisions 
of the projectile p with successive target nucleons j and f is given byi 

((Tpj(Tpj')i = o- 2 (l + u a ) , (47) 

and the scaled variance of the cross-section fluctuations, u a , is given by 

= ^ - I- (48) 

The freezing of the projectile configuration during the passage across the nucleus corre- 
lates <Tpj and a P ji through the internal coordinates of the projectile. When the projectile 
is in a small configuration it has small probability to scatter on both j and j', and when 
in a large configuration is has a large probability to do so. [On the other hand, a slow 
projectile has ample time to change its internal structure between collisions; subsequent 



We assume that the internal degrees of freedom decouple from the nuclear wave function 
coordinates, i.e., we neglect small modifications of the properties of nucleons in nuclei as 



collisions are uncorrelated and then (cr p jCr p j>)i = cr 2 , so that uj a = 0.] As we showed in 
Refs. [j5|,fF yi2t|l3p3| , the size and form of the cross-section fluctuations can be extracted 



from diffractive scattering experiments. 

Fluctuations of the color degrees of freedom of a hadron lead to fluctuations in the 
overall scale as well as shape and range of its interaction cross-section da/d 2 b with other 
hadrons. In the following we shall assume for simplicity that only the scale a of the 
cross-section fluctuates but not its shape or range r a [see, e.g., Eq. (0)]. Fluctuations 
in shape and range can lead to further corrections to the effects we found from including 
the finite interaction range, which should be modelled in the future. Employing the 
parametrization (|24]), we see that a 2 enters as a common prefactor in (0); inclusion of 



cross-section fluctuations at this level thus simply multiplies the right side of ([?[) by a 
factor 1 + u> a , and we directly obtain the scaled variance in proton- nucleus collisions, 

uj n = l-a + oj a ((N) -a); (49) 

the important new ingredient is the last term due to the fluctuations in cross sections. For 
a central p- 208 Pb collision, (N) ~ 7, and with uj a ~ 0.2-0.3 at high energies as estimated 
in Ref. cross-section fluctuations contribute ~ 1.5-2.0 to un, thus increasing by 
a factor 2-3. 

The fluctuations in nuclear collisions can be calculated analogously to the proton- 
nucleus case. In the four terms of Eq. (p9|), only the second and third involve multiple 
collisions of the same nucleon; therefore in these terms we must replace the prefactor a 2 
by (x 2 (l +u> a ). The fluctuations in the number of NN collisions in central nucleus-nucleus 
collisions for £ < Rb are thus given by 

uj n = 1 + 7 - a - (5 + u a {N pA + N pB - a - (3) . (50) 

We can interpret the covterm as the correlation coming from each projectile nucleon 
making on average iVp^ — a hits in the target, and a target nucleon being hit on average 
N P B — P subsequent times by projectile nucleons. Equation (|50|) generalizes the proton- 
nucleus result of Eq. ( f49D to nuclear collisions by adding the correlations from target 
nucleons being hit multiple times. Expressions for N p a and iVp£ are given in the previous 
section, Eqs. P])-©. 



V. TRANSVERSE ENERGY SPECTRA 



nucleus-nucleus collisions use either the number of binary NN-collisions or the number 
of participating nucleons for the sources of particle or E? production [p4}| . We first 
discuss the fluctuations in a binary collision picture, and then turn to the participant, 
or "wounded," nucleon picture. 



A. Fluctuations in Binary collisions 

In an independent-source description of collisions, the £r-distribution is obtained 
by folding the Et contributions from N sources, where the distribution of the number 
of sources has mean (N). The overall form of the transverse energy distribution in a 
high-energy nucleus-nucleus collision is determined by the geometry of the collision, as is 
well-illustrated in the NA34 geometric model, where the Et distribution is obtained by 
folding the Et contribution from N sources (which becomes Gaussianly distributed for 
sufficiently large N), multiplying by the distribution of sources da/dN, and summing 
over N: 

d(T r d(J e -(E T -Ne )y2uNel 

~d~Er~ J dN {2nujNe 2 y/ 2 ' (51) 

The parameter eo is interpreted in the NA34 model as the average Et produced per 
source, while u describes the standard deviation or fluctuation around the mean value; 
both parameters are extracted by fitting to data. We do not take this interpretation 
literally here, but rather adopt the point of view that the model provides a useful 
parametrization of the data. 

The cross section for making N sources, which is taken to be the number of binary 
collisions, is 

^ = J d 2 bP N (b), (52) 

where -P/v(fe) is the probability for making N sources in a collision at impact parameter 
b. In the simplest geometric model, P/y is non-zero only at the mean number of sources, 
N(b), for given impact parameter; for example, for small projectiles B and spherical 
target nuclei of radius Ra, one has N(b) = 2Bap (R 2 A - b 2 ) ' . More generally one can 
take P/v(6) to be a Poisson distribution around the mean, but, as we have discussed 
above, correlations lead in fact to a much narrower structure for P^ than Poisson. 
Equation (|5l"l) also describes the multiplicity distribution when E T is replaced by the 
multiplicity n; the parameter eo is then replaced by the average number of particles 



fluctuations. The effect of impact parameter averaging is in fact separated out in the 
NA34 analyses of nucleus-nucleus collisions f| . As described in Ref. |25j , the tail of the 



Et distribution in ultrarelativistic heavy-ion collisions is determined by central events, 
on which we concentrate here.! The NA34 analysis thus provides a reliable, albeit 
model-dependent, way to extract the fluctuations uj for central (b = 0) collisions from 
data, and yields values of to in collisions of 32 S on nuclear targets which range, in the 
pseudorapidity interval —0.1 < rj < 2.9, from ~ 1.5 to ~ 5.8. The uj's extracted from the 
NA34 analyses directly describe the transverse-energy fluctuations in central collisions, 

_ /A rt ( E t) ~ ( E t) 2 /„\ 

u=(N) , (53) 

where (N) is evaluated at zero impact parameter. The problem is to understand the 
physics producing the large values of uj needed by NA34. As we first discussed in Ref. 
H, and expand on here, cross-section fluctuations lead to a substantial contribution to 
oj, capable of explaining deviations of standard estimates from experiment ||||. 
The principal contributions to uj are: 

uj = uj + uj def + uj n , (54) 

where uj is the width of the single source .^-distribution, uj de f = 4(iV)<5 2 /45 is the 
variance of the fluctuations from the deformation of the target, with 5 the nuclear 
deformation parameter, particularly large for W and U. The third term, ujn = {(N 2 ) — 
(N) 2 ) I (N), takes into account fluctuations in the number of sources in a central collision. 
With the assumption that the number of sources is given by the number of binary 
collisions, fluctuations in the number of sources is given by Eq. (|35f) in the absence of 
cross-section fluctuations, and Eq. (|50"D with cross-section fluctuations included.0 



The procedure of extracting fluctuations at given impact parameter is much more difficult 
for proton-nucleus collisions which do not exhibit the plateau spectrum in nucleus-nucleus col- 
lisions. While the characteristic "shoulder" in nucleus-nucleus -Er-spectra marks the beginning 
of region dominated by central collisions, in proton-nucleus collisions impact-parameter aver- 
aging results in a long tail which is further smeared by fluctuations as well as degradation and 
rescattering p|. 



Using u m 0.5 as determined from pp-collisions in the comparable energy and rapid- 
ity range, one cannot, with these assumptions in the absence of cross-section fluctuations, 
explain the experimentally observed large fluctuations.0 With a ~ 32 mb and po — 0.15 
fm -3 , and a heavy target we have N p a ~ 2ap R,A ~ 1.2A 1 ^ . In CERN and Brookhaven 
experiments, with 16 or 32 S projectiles, N p b — §(XPo-Rb — 2 — 3. We show in Fig. 4 
uj — uj de f = un + ^o given by Eq. (|50"D , for 32 S, with various values of u a , taking uj q = 0.5, 
a — (3 — 0.25 and 7 = 0.15, as discussed above. We see that cross-section fluctuations 
are able to account for the large a>'s found experimentally with u a ~ 0.25, a value con- 
sistent with that extracted from inelastic shadowing in nucleon-deuteron scattering as 
well as forward diffractive scattering amplitudes at CERN-SPS energies [fL2"] . Equation 
(^0|) describes the 16 O data |28| as well with similar values for uj a . To disentangle the 
color transparency and opacity effects from other sources of Et fluctuations in future 
experiments, especially rescattering, it will be useful to study the energy and rapidity 
dependence of the E T fluctuations. 



B. Fluctuations in Participants or "Wounded nucleons" 

Let us consider briefly the fluctuations in a model in which the number of sources is 
the number of participant or wounded nucleons. The number of sources in proton-nucleus 



of GeV), which leads to uj res = (0.45) 2 (iV) /Ec e v of order 0.2; we are grateful to Dr. H. 
Strobele for emphasizing the existence of this contribution. In Eq. (|54|) we also have neglected 
terms due to decrease of the transverse-energy production per source in successive hits. As 
argued in Ref. || these terms lead to a small reduction of oj. 

11 Rescattering of secondaries might be expected to broaden the tails of the -E/r-distribution 
and thus increase u, since it allows a secondary to lead to a wider range of transverse energy and 
multiplicity in final states. But rescattering also increases the number of scatterings and thus 
tends to narrow the fluctuations; in general, adding more degrees of freedom to a system while 
keeping its average output fixed decreases its fluctuations. A calculation for central collisions 



of 32 S on Au at CERN-SPS energies by Werner fl27| using the VENUS 4.16 code both with 



collisions is very similar for the binary-collision and the wounded-nucleon pictures, since 
the number of participants is just the number of binary collisions plus one, the projectile. 
Similarly, the fluctuation around the mean number of sources is almost the same. For 
nuclear collisions the numbers of sources in the two pictures can differ considerably fll! 



The geometric-model result (pip can also be applied in the wounded-nucleon picture. 
The interpretation of the source parameters eo and ujq is, however, very different in the 
two pictures. In the binary collision model eo is the transverse energy produced in a 
nucleon-nucleon collision whereas in the wounded nucleon model it is half that value. 



As we see from Eq. (|53|), f2 = u/{N) is the same in both cases. 



The number of participants in a nucleus-nucleus collision is given by 

(N) WN ~ J d 2 s [T fl (s)(l - e - CTT ^ s+b >) + T A (s + b)(l - e"^)] , (55) 

a formula which assumes a <C -nR 2 ^. When the projectile is much smaller than the target 
nucleus this result simplifies to 

(N) WN ~ B{\ - e- N ^) + irR 2 B T A (b) ^1 - J^[l - (1 + N pB )e~^] \ , (56) 

where N p a = crT^(6) and N pB = aT B (0). While N p a is the same as in Eq. (|38"D, N pB is 
a factor 3/4 smaller than N pB in Eq. (|39]) . 

In the limit in which all nucleons in the overlapping volume of the two nuclei in the 
collision participate, i.e., small nucleon mean free-paths, corresponding to large N p a and 
N p b, and assuming sharp nuclear surfaces we find from (|55|) that for central collisions 



(N) WN = B + A-\A 2 I*-B 2 '*f\ (57) 

For small projectiles (|5T| ) reduces to 

(AO WN ~ B + 7T R\T A (b) , B < A , (58) 

valid in this form also for non-central impact parameters. In comparison, the mean 
number of binary collisions, (A) ~ BaTA{b), is generally is larger. 

The latter results are, to first approximation, independent of the cross section; there- 
fore we do not expect fluctuations in cry in this limit to affect the distribution of wounded 
nucleons significantly. 

To estimate the fluctuations in the number of participants including correlations, we 
limit ourselves for simplicity to central collisions and assume sharp nuclear surfaces and 



(N(N - l)) WNeA = J d 3 rd 3 r' p%\ r ,r')6(R B - \r ± \)6(R B - |r' ± |) . (59) 

[Note the formal similarity of this equation to Eq. (§) for pA collisions with the range 
r s of the nucleon-nucleon interaction replaced by Rb] For an uncorrected nuclear 
gas, for which £ = Ra and C 2 = 1/A, Eq. (|59"D gives = (N)^r NeA (l — 1/A) and so 
ojn = 1 — (N)/A, the usual binomial result (cf. |14[)) . More realistically, £ <C i n which 
case Eq. (|59|) implies (N 2 )wNeA = (^)wNeA an d thus cjjv = (cf. (|18D ). Correlations 
effectively eliminate fluctuations in Nwn- Fluctuations in observables arise only from 
fluctuations in the individual sources, given by ujq^wn = 2^0, since two wounded nucleons 
make up one binary NN collision. To compare the predictions of fluctuations in the 
wounded-nucleon picture to the NA34 data, we should keep in mind the conventional 
but arbitrary scaling by (N) in the relation of fl53|) to the fluctuations in transverse 
energy, and thus compare ujwn = ^o,wn{N) / (N)wn — 1 — 2 to the u extracted in the 
NA34 geometric model. The fluctuations in the wounded-nucleon picture are generally 
not large enough to explain the NA34 data. 

VI. CONCLUSION 

As we have seen, nucleon correlations in projectile and target nuclei, the range of 
the underlying NN interactions, and cross-section fluctuations are all important for de- 
termining fluctuations in the number of NN collisions and thus observables in nuclear 
collisions. Correlations reduce the density fluctuations and therefore also the fluctu- 
ations in the number of NN collisions. Because the nucleus is a self-bound liquid of 
nucleons, correlations are particularly important in nucleus-nucleus collisions, and lead 
to significant reduction in fluctuations compared to those for a gas of nucleons assumed 
freely swarming in the nucleus. Interestingly, the correlations in the incident nuclei are 
one of the few features of the original nuclear structure that are a determining factor 
in high-energy collisions. Allowing for a larger range for the interaction or collision 
probability reduces the effect of the correlations between two nearby nucleons, and in- 
creases the fluctuations. Cross-section fluctuations are an important feature of collisions 
of a hadronic projectile at sufficiently high energy that time dilation does not allow the 
projectile to change its configuration between successive collisions in the target nucleus. 
In this case successive collisions become correlated, with the effect of increasing the 
fluctuations in the number of NN collisions significantly. 

We emphasize the need to incorporate the physics we discussed here into Monte-Carlo 



separated by a minimum distance does not reproduce nuclear correlations and fluctua- 
tions satisfactorily, especially in nucleus-nucleus collisions; approaches to implementing 
nuclear correlations that correctly describe a nuclear liquid were discussed at the end of 
Sec. 3. To a first approximation, cross-section fluctuations can be incorporated in the 
probabilistic language of Sec. 3 by assuming factorization for interactions of two nu- 
cleons (1 and 2) in different configurations, i.e., assuming the cross section to have the 
form, a = a\Oija (cf. discussion in |||7j]), and then generating a for each of the nucleons 
with individual weights P{<7\) and P(cr 2 ), where P(u) is the normalized distribution of 
cross sections for NN scattering fll2"f . 

Beyond fluctuations in the probability of collisions, one should also incorporate fluc- 
tuations in the outcomes of individual collisions, i.e., fluctuations on the parton level. 
Since parton distributions depend on the size of the nucleon configuration, fluctuations 
in internal nucleon configurations lead directly to fluctuations in parton distributions 
and hard parton-parton interactions. While microscopic collision models in terms of 
partons have been developed for ultrarelativistic heavy-ion collisions (see, e.g., [^ ,[30H ) , 
they so far only take fluctuations in number of collision partners into account, employing 
average structure functions, or parton momentum distribution functions, for individual 
nucleons. A simple way to take fluctuations on the parton level in account in such mod- 
els, based on the fact that increase of the effective transverse size of a hadron allows 
emission of softer partons, is to let the parton distributions Pn{x, Q 2 , c) scale as 

p N (x,Q 2 ,a) = p N (x,Q 2 a/a), (60) 



similar to the rescaling model of the EMC effect 

At small x one also has to take into account shadowing and enhancement effects which 



depend both on a and the impact parameter of the colliding nucleon [ 32 1 ; discussion of 
these effects is beyond the scope of this paper. Furthermore, observed fluctuations in 
hard parton-parton interactions will depend on the subset of events on which one triggers 



in a collision (see |33| and discussion in e.g., biasing towards high events will 

enhance the probability of large interaction cross-section configurations occurring in the 
colliding nuclei. Methods to study this phenomenon experimentally in pA collisions are 



discussed in [34 . 
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APPENDIX A: THE AGK CUTTING RULES AND BINARY COLLISIONS 

The AGK diagram cutting rules for high-energy scattering processes [||,|| lead to the 
theorem that in nucleon-nucleus (NA) scattering, inclusive particle production away from 
the projectile fragmentation region is simply that in NN scattering times the number 
of nucleons, A, in the target nucleus. In this Appendix we briefly discuss the extension 
of this result to nucleus-nucleus collisions at central rapidities, to give an indication 
of how the probabilistic methods we have used in this paper can be extracted from a 
more formal field-theoretic approach. The prediction of the theorem for nucleus-nucleus 
collisions that particle production scales with the number of NN collisions supports the 
use in this paper of the binary collision model in describing fluctuations in the central 
rapidity range. 

The theorem is based on the Glauber multiple NN scattering formalism within the 
eikonal approximation, and thus ignores rescattering of produced particles during the 
collision. Rescattering of target nucleons on other target nucleons is not taken into 
account, and in nucleus-nucleus collisions neither is scattering of projectile nucleons on 
other projectile nucleons. In addition, coherent processes, such as showers of particles 
overlapping spatially and interacting strongly, are ignored, cLS IS cL possible transition 
to a quark-gluon plasma in which particle production cannot be described in terms of 
nucleon-nucleon interactions. Finally, the AGK theorem ignores constraints from global 
energy conservation in multiple NN collisions, and is thus valid only at those rapidities 
and at sufficiently large energies that particle production depends slowly on the initial 
energy. Within these assumptions, AGK show that there is an exact cancellation of 
processes involving inelastic interactions of the projectile with more than one nucleon at 
a time. For example, in the interaction of the projectile with two nucleons, the forward 
amplitudes for the two processes where the projectile interacts inelastically with one of 
the nucleons and elastically with the other cancels against the amplitude for production 
of particles in inelastic interactions with both nucleons. The net result is that the total 
multiplicity in inelastic collisions is the sum of the multiplicities in inelastic collisions of 



constrains particle production; the cancellation fails in the projectile fragmentation re- 
gion because the way in which the projectile energy is shared in inelastic collisions 
depends on the number of target nucleons with which it collides. The net result is 
shadowing of particle production close to the projectile rapidity.S The more nucleons 
involved in a collision, the stronger is the constraint of global energy conservation, and 
the wider the rapidity interval where the AGK cancellation is not effective. For exam- 
ple, in central nucleus- nucleus collisions at lab energy of 200 GeV/A, as at CERN, the 
limited total energy constrains particle production in the nuclear fragmentation regions, 
and precludes application of the AGK theorem there. 

The extension of the AGK theorem to nucleus-nucleus (BA) collisions, for particle 
production away from nuclear fragmentation regions, is: 

a BA ^ a+x = BAa NN ^ a+x , (Al) 

where ( j BA - >a + x an d a NN ^ a + x are the inclusive cross sections for producing particle a 
in BA and NN collisions respectively. The number of particles a produced per collision 
is 

BA-^a+X 

- BA ^ = —jnr- > ( A2 ) 

where a BA is the total production cross section for colliding nuclei A and B (i.e., not 
including quasielastic processes in which no hadrons are produced); for NN scattering 
a NN = a. The mean number of NN collisions averaged over impact parameter is then 

BA-^a NN 
( N ) = ^NN^a = BA ^BA> (A3) 

a result in fact independent of the particle a. By writing this equation as 

( N ) = J^A ° NN j d2sT ^ s + h ) T B(s), (A4) 



we may interpret Eq. (|A3|) as an average of impact parameter of the number of collisions 
at given impact parameter, 

(N(b)} = a NN J d 2 sT A (s + b)T B (s); (A5) 



One can, in fact, extend the AGK arguments to include energy conservation within the 



cf. Eq. (p7|) . The generalization of Eq. (|A4j) to a finite interaction range is outlined and 
applied to diffractive scattering and inelastic shadowing effects in Refs. p6|j7|Jl^| . 

The above example indicates the correspondence between the probabilistic language 
of Sec. 3 and field theory. The AGK theorem applied to nucleus-nucleus collisions leads 
to the same partial probabilities as the binary collision model of collisions involving a 
given number of nucleons in inelastic interactions. Binary collision models make par- 
ticular assumptions about how energy is shared in the various sub collisions. At central 
rapidities, where the energy constraints are unimportant, the results are independent 
of these assumptions, and one can apply the AGK theorem. However, binary collision 
models give specific predictions in the fragmentation regions as well. 

At very high energies the contribution of minijets to particle and production 
starts to become comparable to that from soft processes. According to Eskola et al. 
PH minijets contribute more than 50% at energies above RHIC energies, A/i^200GeV 
(such estimates, however, carry the uncertainty of how low in transverse momentum one 
can use perturbative qcd to calculate jet production). Since minijets are hard processes 
with kinematics where nuclear shadowing in the parton distributions can be neglected, 
they are produced independently in each NN collision. Hence in this case the binary 
collision model for parton-parton scattering should be valid, although one has to take 
into account fluctuations in the parton distributions on nucleons, cf. Eq. fl60"|). 

The inclusion of color fluctuation effects in the description of nucleus-nucleus colli- 
sions for calculation of total cross sections is equivalent in qcd to inclusion of Gribov 
inelastic shadowing effects, the presence of inelastic intermediate states in the inelas- 
tic eikonal approximation (at the level where one neglects the high-mass diffraction, 
described by triple Pomeron diagrams). The use of the eigenstate formalism |37| , |3"8|j (ap- 
plied in this context in Refs. |I2|,23]|) allows one to diagonalize the transition matrices in 
the inelastic eikonal diagrams and write these diagrams as a sum of terms corresponding 
to the scattering of nucleons in the states with fixed a. The probabilistic interpretation 
of the eikonal approximation for the multiple scattering combinatorics is thus preserved 
even in the presence of color fluctuations, allowing one to understand the results of Sec. 
3 from a field-theoretic perspective. For example, Eq. (29), with color fluctuations, 
can be derived using the AGK rules for the double inclusive spectrum. Indeed, any 
quantity calculated in the eikonal approximation proportional to a n is given, with color 
fluctuations taken into account, by the same expression with the substitution o n — > (<J n ). 

We also note that the AGK cutting rules connect the cross sections for inelastic in- 
coherent hadron production and diffractive hadron production, as arising from different 



inelastic incoherent hadron production. The agreement of the results of the calculation 
of inelastic coherent diffractive processes in hadron-nucleus scattering with existing data 
4C)| provides a further confirmation of the picture of color fluctuations adopted here. 
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FIGURES 

FIG. 1. Straight-line geometry of a proton- nucleus collision displaying the length scales of 
the cross-sectional distance, r s , the interparticle spacing, tq, and the correlation length, £. 

FIG. 2. Correlation functions from j[{J (full curve), the Gaussian fit, Eq. (|^) (dashed 
curve), an ideal gas, Ca = 1/A (dotted curve) and a gas with excluded volume (long-dashed 
curve) for a nucleus with A = 208. 

FIG. 3. Fluctuations ujn for central collisions of 32 S on different targets for a strongly 
correlated nucleus, Eq. ( |50j ) (solid line), as well as for a nuclear gas, Eq. (|37D (dashed line). 

FIG. 4. Fluctuations uj — LOdef = ujq + lun for central collisions of 32 S on different targets 



for uj a = 0.0, 0.1, 0.2 and 0.3 calculated with Eq. (50) . The experimental values (dots) are 
taken from Ref. 01. 



